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DECODING SYSTEM AND METHOD FOR DIGITAL COMMUNICATIONS 

This is a continuation of part of U.S. Application No. 09/535,902. Priority is claimed 
based on said Application No. 09/535,902 and based also on Provisional Application No. 
5 60/271,615 filed February 26, 2001. 

Background Of The Invention 

1 . Field of the Invention 

The present invention relates to digital communication systems. More particularly, the 
invention relates to decision decoding of received digital information for use in subsequent 
10 decoding processes as well as to weighting selected received coded data based on its 
estimated reliability prior to decoding. 

2. Background of the Related Art 

5 The process of transmitting digital information may in one sense be thought of as a four- 

? step process as shown in FIG. 1. First, an incoming data stream of electrical digital signals 
11 5 x(t) is coded by a coder 1 0 to produce coded information C(x(t)) that is provided to a trans- 
it mitter 20 which transmits the coded information C(x(t)) over a communication medium 30. 
3 The coder 10 processes the data stream x(t) so that the coded information C(x(t)) is robust to 
3 errors; that is, errors in the coded information C(x(t)) can be detected and corrected. This is 
H necessary because, as will be seen, the transmission process introduces a noise component 
y20 into the transmitted signal which would corrupt the signal and render it useless were it 
J impossible to be removed. 

In an exemplary case, the incoming data x(t) is digitized audio such as voice or music 
from a known source; the coding unit 10 applies a coding transformation to the incoming data 
x(t) to produce the coded information C(x(t)); and the transmitter 20 modulates a radio- 
25 frequency signal with the coded data C(x(t)) to send it as a radio-frequency signal through the 
atmosphere which serves as the communication medium 30. 

In traversing the communication medium 30, the transmitted information C(x(t)) 
acquires a noise component n(t), and the resultant signal C(x(t)) + n(t) received by a receiver 
40 is decoded by a decoding unit 50 to recover the original data stream. Through the use of 
30 error-checking and correction properties of the coding algorithm employed by the decoder 
50, the noise component n(t) can be detected and removed from the received signal 
C(x(t))+n(t), after which the signal C(x(t)) is processed in the decoder 50 using an algorithm 
complementary to the one used in the coder 10 to obtain the original data stream x(t). 
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The term "noise" as used herein is intended to include any source of signal degradation, 
including but not limited to thermal noise, circuit noise and interference. For single carrier 
(equalized) systems, the preferred scheme includes identifying those bits in particular bit 
patterns that are most likely to be received in error, and weighting/puncturing those bits. 
5 In the above example, the receiver 40 demodulates the received radio frequency signal to 

produce the received coded information C(x(t))+n(t). The received coded information is 
error-checked and decoded by the decoding unit 50 to produce the outgoing digital data 
stream which might then be converted to an analog signal to be amplified for driving a 
speaker. 

10 It is common in systems such as the one described above for the coder 10 to code C(x(t)) 

using a convolutional coding system* In contrast to a block code where a block of bits is 
coded according to a rule defined by the code without regard to previous inputs to the coder, 
;|> a convolutional code additionally makes use of previous bits it processed in the coding 
; ^ process; thus, the coder 10 is a type of finite state machine. 

(15 An example of a convolutional coder of a type which might be used in the coder 10 is 

;U shown in FIG. 2A, As is evident from the coder block diagram, coding in this unit depends 
^ 3 not only on the state of the current input bit x n but also on the state of the two previous input 
i 5 bits held in the two latches D. Consequently, the coder is a four-state machine as shown by 
I .] its state diagram in FIG. 2B. From the state diagram it is relatively straightforward to derive 
I L ^20 the coder's trellis diagram as shown in FIG. 2C. Here, the branch to a particular state 

corresponding to a zero input is represented by a solid arrow, and the branch to a particular 
state corresponding to a one input is represented by a dotted arrow. 

With the trellis diagram of FIG. 2C, the coding process will be readily apparent to those 
skilled in the art. Applying an input vector of, e.g., X = (1101011) produces outputs from 
25 the coder of FIG. 2A as shown in TABLE I (assuming the coder starts from an all-zero state). 



60215797vl 



-3- 



ATH-021 cip 



Input 


Internal State D X D 2 


Output yiy 2 


Xl 










1 


0 


0 


1 


1 


1 


1 


0 


0 


1 


0 


1 


1 


0 


1 


1 


0 


1 


0 


0 


0 


1 


0 


1 


0 


1 


0 


1 


0 


0 


1 


1 


0 


0 


1 


0 


1 


1 


0 


1 


0 


0 


1 


1 


1 



TABLE I 

Thus, for the coder of FIG. 2A, C(X) = (1 1101000010010101 1). 
"5 s 5 If the communication medium 30 were error-free and the receiver 40 were assured of 

; & providing C(X) to the decoder 50, it could recover the original data stream simply by 
r ;f; applying C(X) to the trellis network of the coder state machine and track the path 

therethrough which generated the sequence C(X). For example, C(X) = (1 1101000010 
? 0 10 10 1 1) would be generated by the path shown in FIG. 3, where states occupied by the 
; 3IO coder are shaded and branches taken by the coder are in bold, with the corresponding coder 
output C(X) given above each stage. From this, the original input sequence X = (l 10101 
•M 1 ) can be obtained. 

'Cl The above discussion and those hereinafter assume that the coder 10 and decoder 50 

operate according to the same coding algorithm; that is, they both base their processing on the 

15 same coder circuit, state machine and the like; thus, dimensions of the trellis network 

necessary to accurately represent the coding algorithm are known to the decoder, the number 
of stages necessary in the trellis network is known, etc. Further, it is assumed that the coder 
state machine begins in an all-zero state, and that the input stream X is padded with a 
sufficient number of zeroes to return it to such a state at the conclusion of the 

20 coding/decoding process. 

Assume now that the signal received and processed by the receiver 40 has a non-zero 
noise component n(t), i.e., C'(X) = C(X) + n(t). This might result in the decoding unit 50 
receiving C'(X) = (11101700017010101 1), with the sixth and eleventh bits being 
errors due to n(t). In this case, the original input sequence cannot be found simply by 

25 applying C'(X) to the trellis of FIG. 2C. A typical error-correction algorithm used in this 
situation is the Viterbi error correction algorithm described as follows. 
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1 . Given the trellis network, associate a metric with each stage in the trellis and set the 
metrics for all states in the first stage to zero. 

2. For each state in the next stage, find the "distance" of each branch to it from the next- 
received subsequence of C'(X) and add it to the metric of the current state from which it 

5 branches. 

3. Choose the minimum of the values calculated in Step 2 as the metric for the given 
stage in the next stage, and choose the branch leading to it as the survivor branch to that state. 

4. Repeat Steps 2 and 3 until the end of the trellis network is reached. 

5. Select the state in the last stage of the trellis network having the minimum metric and 
10 work backwards to the beginning by selecting survivor paths to produce the best guess about 

C(X). 

; r If the coder is known to start from an all-zero state, the system is preferably designed so 

that traceback is always done to the all-zero state in the first stage. This can be done by, for 
example, initializing the all-zero state to zero and initializing the other states to large values. 

, J. 1 5 Also, there are various techniques available to determine the state from which to begin a 

f;* traceback. For example, at the end of a frame when the coder ends in a known state (in the 
example, 00), start the traceback from that known state. 

: t Application of the Viterbi decoding algorithm to the corrupted stream C'(X) and the 

trellis network of FIG. 2C is shown in FIG. 4. Here, 

; %20 — the metric for a state is shown at its center; 

— the distance of a current subsequence from a given branch is shown by a number 
above the branch; 

~ survivor branches (except for those in the best guess path) have a white arrowhead; 

and 

25 — the best guess path is in bold. 

First, metrics of all states in the initial stage are set to zero. Then, the distance (here, the 
Hamming distance) between the subsequence 1 1 and each branch value is assigned to that 
branch. For the first node, the two branches are 00 and 1 1, so they receive values of 2 and 0. 
The next node has branches 01 and 10, so each of these branches receives a metric of 1 . The 

30 third node has branches 1 1 and 00, so they respectively receive values of 0 and 2. Finally, 
the fourth node has branches 10 and 01 which each receive a Hamming distance of 1. 

Then, for the first state in the second stage, the two branch values to it are 2 and 0, so its 
metric is set at 0 and the 0 branch from the third state in the previous stage is set as its 
survivor branch. For the second state in the second stage, the two branch values are 0 and 2, 
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so its metric is set at 0 and the 0 branch from the first stage in the first stage is set as its 
survivor branch. Next, the third stage in the second stage receives values of 1 and 1, so its 
metric is set to 1 and one side is randomly or deterministically set as the survivor branch. 
Finally, the last state in the second stage receives values of 1 and 1, so it has a metric of 1 and 
5 both branches are set as survivor branches. 

This process is repeated to the end of the trellis as shown in FIG. 4. Then, starting at the 
first node in the last stage (since it has the lowest metric of the last stage), survivor paths are 
traced back to the first stage. As can be seen from FIG. 4, the result is the same path traveled 
in the error-free case. Thus, the Viterbi algorithm has successfully eliminated the error 
1 0 contribution of n(t). 

Although the Viterbi algorithm works well for situations such as those given above, it is 
f ^ not always practical to directly implement the algorithm in this way. For example, some 
>f communication systems make use of alternative information-carrying techniques such as 

quadrature amplitude modulation (QAM), which requires alternative methods of retrieving 
:;|Ji5 the received data stream. Additionally, some communication systems make use of 

interleaving and other techniques which restricts the type of decoding systems that can be 
; / used therein. Thus, there is a need for a decoding technique that addresses these issues. 
* Also, although the Viterbi algorithm works well for the situations described above, there 

■^""',■3," 

■ L r ; are improvements that can be made to increase its usefulness. For example, consider the case 
^20 of an OFDM system with interleaved data mapped to multilevel QAM constellation points. 
^ In this case, it is impossible to use Euclidean distance measurements to determine the proper 
weights to use in the algorithm, since the multiple bits represented by a single measurement 
are non-sequential in the trellis. 

Further, the Viterbi algorithm in its current state does not allow for unequal weighting of 
25 bits based on knowledge of their error probability. This unequal error probability may be due 
to relative position in a multi-level constellation, or position in time, frequency, or space 
which is more or less likely to be susceptible to noise or interference than other positions. 

Also, one of the primary burdens for wireless data transceivers, in terms of both 
processing latency and consumed integrated circuit real estate, is data decoding such as 
30 described above. Available wireless receivers fail to differentiate between reliable (good) 
and unreliable (bad) bits. Thus, even though bad bits are more difficult to decode and are 
often discarded once decoded, all bits, both good and bad, are presently utilized equally in 
available decoding processes. There is therefore a need in the art for a method of analyzing 
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data to determine, in real time, its reliability and utilizing such information before decoding 
the analyzed data. 

Summary Of The Invention 
In view of the above deficiencies of the prior art, it is an object of the invention to 
5 provide a digital information decoding system which can be used in a wide variety of 
communication systems. 

It is another object of the invention to provide a digital information decoding system that 
can be used with a wide variety of information demodulation techniques. 

It is a further object of the invention to provide a digital information decoding system 
10 that can be used with interleaving digital communication systems. 

It is yet another object of the invention to provide a technique for analyzing data to 
Z' determine, in real time, reliability of received data and using such information in decoding 
}f: the analyzed data. 

^ The above objects are achieved according to a first aspect of the present invention by 

;S1 5 providing a Viterbi decoding system which interprets bits in received QAM constellations as 
many-valued parameters rather than binary valued parameters. It performs the Viterbi 
algorithm using these many- valued parameters to provide results superior to hard decision 
y decoding. Rather than applying a hard 0-1 function to the QAM data, the system uses a non- 
b 5 stepped linear or curved transfer function to assign values to the bits. 

'120 According to another aspect of the invention, these objects are achieved by providing a 

system which differentiates between data bits based on their estimated reliability, giving 
more emphasis to decoding reliable bits than unreliable bits using any of a variety of 
techniques. By differentiating between good and bad bits and de-emphasizing or ignoring 
unreliable bits, the system can provide a significant reduction in uncorrectable errors and 
25 packet loss. 

According to another aspect of the invention, this process of selective weighting or 
"puncturing" enables the discarding ("hard puncturing") or de-emphasizing ("soft 
puncturing") of received bits deemed less reliable than a particular minimum threshold level. 
Within a multi-carrier environment, such a scheme further involves identifying which carriers 
30 have poor signal-to-noise ("SNR") characteristics and attendant high error rates, and 
puncturing/weighting bits in the identified carriers. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

These and other objects, features and advantages of the present invention are better 
understood by reading the following detailed description of the preferred embodiment, taken 
in conjunction with the accompanying drawings, in which: 
5 FIG. 1 is a block diagram of a digital communication system according to the prior art; 

FIG. 2A is a diagram of a convolutional coder used in the system of FIG. 1; FIG. 2B is a 
state diagram of the coder; and FIG. 2C is a trellis network of the coder; 

FIG. 3 shows use of the trellis network in decoding an error-free output data stream from 
the coder of FIG. 2; 

10 FIG. 4 shows use of the trellis network in decoding an output data stream from the coder 

of FIG. 2 having errors; 

FIG. 5 shows a QAM constellation according to the prior art; 
3 FIGS. 6A - 6D show bitwise decision decoding functions for the QAM constellation of 

I FIG. 5; 

jll5 FIGS. 7 A and 7B show transfer curves corresponding to the decoding functions of FIGS. 

6A and 6B; 

FIGS. 8A and 8B show transfer curves according to a first preferred embodiment of the 
^: present invention; 

■ U FIG. 9 shows use of a trellis network in decoding a QAM output data stream having 

:il0 errors; 

^ FIGS. 10A and 10B show transfer curves according to a first preferred embodiment of 

the present invention; 

FIG. 1 1 shows channel frequency response with blocker signals or other interference 
present; 

25 FIG. 12 is a block diagram of a transmitter in an OFDM system; 

FIG. 13 is a block diagram of a receiver in an OFDM system; 

FIG. 14 shows a portion of the Viterbi decoder according to the present invention; and 
FIGS. 15A and 15B illustrate a decision directed feedback according to the present 
invention. 

30 Detailed Description Of Presently Preferred Exemplary Embodiments 

The above arrangement described in the Background section works relatively well if the 

digital information is used to modulate the carrier signal directly; however, this is not always 

the case. For a variety of reasons, some communication systems employ alternative arrange- 



i 
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ments. For example, the IEEE 802.1 la communication standard uses a quadrature amplitude 
modulation (QAM) system in which the modulated carrier signal does not directly correspond 
to the information being transmitted, but to a position in a two-dimensional array or "constel- 
lation" of possible binary values. This constellation is known to the sender and receiver, so 
5 that a position transmitted to a receiving station can be accurately correlated with a corres- 
ponding digital sequence in the constellation. 

Consider FIG. 5. This is a constellation which might be used in a 16QAM (16 bit 
quadrature amplitude modulation) coding system. Here, rather than providing a single digital 
output to the decoder 50, the receiver 40 would generate two signals I and Q. So, a data point 
10 (I, Q) P of, e.g., (1 , -3) would correspond to the binary sequence 1 1 00. This sequence 1 1 00 
could then be applied to the Viterbi trellis as a received subsequence as described above. 

In this case, noise effects from the communication medium 30 could vary I and Q so that 
C they do not point directly to the appropriate constellation point. For example, rather than (1, - 
:I 3) the receiver 40 may generate a data point P 5 of (1 .25, -2.75). This may be dealt with using 
; vl 5 either a hard decision decoding process or a soft decision decoding process. 
^ In the soft decision decoding process, a Euclidean distance measurement is used to 

determine the constellation point corresponding to the received (I, Q) point. For P' = (1 .25, - 
f 2.75), the distances are 

J D noo = ^{1.25 -\f + (-2.75 + 3) 2 * 0.35 

1 Aooo = yJ(l.25-3) 2 +(-2J5 + 3) 2 * 1.77 

D im = yj(1.25-l) 2 ^(~2J5 + l) 2 * 1.77 
Aooi = V(l-25-3) 2 +(-2.75 + l) 2 « 2.47 
20 and so on with increasing distances. This distance would then be used in branch metric 
computation in the rest of the process of Viterbi decoding. 

In the hard decision decoding process, a separate determination is made for each bit as to 
whether it is a 0 or a L For example, given the 16QAM constellation of FIG. 5, a hard 
decision metric for bO would be I > 0 — > bo = 1 as shown in FIG. 6A, wherein the shaded 
25 portion of the constellation indicates the area in which bit bO should be decided to be a binary 
1. A complete set of hard decision metrics would be 
I > 0 b 0 = 1 (FIG. 6 A) 

|I|<2->bi = l (FIG.6B) 
Q>0-*b 2 = l (FIG.6C) 
30 |Q|<2^b 3 = l (FIG. 6D) 
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The transfer curve corresponding to, e.g., FIG. 6A (I > 0 — ► bo = 1) is shown in FIG. 7 A, 
and the transfer curve corresponding to FIG. 6B (| 1 1 < 2 — > bi = 1) is shown in FIG. 7B. It 
should be noted that in FIGs. 7 A and 7B the curves are shown offset from the graph axes for 
illustration purposes only, and it will be understood that, e.g., FIG. 7A has a step precisely on 
5 the abscissa. It will be apparent that similar curves can be produced for FIGS. 6C and 6D. 

As an alternative to such hard decision decoding, a preferred embodiment of the present 
invention uses linear functions in place of the step functions shown in FIGS. 7 A and 7B. 
Two such functions are shown in FIGS. 8A and 8B. The curve shown in FIG. 8 A is the 
counterpart of the curve for b 0 shown in FIG. 7 A, and the curve shown in FIG. 8B is the 
1 0 counterpart of the curve for bi shown in FIG. 7B. Similar curves apply to b2 and b3. As 
before, the curves have been offset from their respective axes for ease of explanation and 
i; preferably are not so offset in actual use. Also, the angles of the transitional portions (or, 

alternatively, the positions at which the transitional portions begin and end) have been chosen 
arbitrarily for purposes of illustration only, and are preferably selected as will be described in 



'fl 5 greater detail below. 

; A complete set of hard decision metrics can be derived from these graphs thusly: 

I<-l^b 0 = O Q<-l-+b 2 = 0 

I -1< K 1 -» b 0 = (I + 1) / 2 -1< Q < 1 -> b 2 = (Q + 1) / 2 

is? 

=j I > 1 b 0 = 1 Q > 1 -» b 2 = 1 

120 

1 5S -3 -> b, = 0 Q < -3 -+ b 3 = 0 

-3<K-l-*bi = (I + 3)/2 -3<Q<-1 -^b 3 = (Q + 3)/2 

-1<I<1 ->-b, = l -l<Q<l->b 3 =l 

KK3^b,=(3-I)/2 KQ<3^b 3 = (3-Q)/2 

25 I>3-»bi = 0 Q^3^b 3 = 0 

The above formulae may be used as is to implement the preferred embodiment; 
alternatively, the I and Q values could be compared to a lookup table to find the appropriate 
bit value (this approach becomes increasingly preferred as the complexity of the formulae 
30 increases). For example, a single table might be used for bO and b2 and when given an I or Q 
value of, say, 0.75, would indicate that the corresponding bit value is 0.875. 

A complete set of bitwise soft decision metrics could also be derived from these graphs: 
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16-OAM 

INPUT MSB 

K-2 0 

-2 < I< 2 (I + 2)/4 

5 2<I 1 



INPUT LSB 

4 < JXJ 0 

|Ij < 4 1 - il| * 0.25 



64-OAM 

INPUT MSB INPUT MID. SIG. B 

K-2 0 6<|I| 0 

10 -2<I<2 (I + 2)/4 |I|<2 1 

2<I 1 2<|I| <6 1 - (|I|-2)*0.25 



g INPUT 

■'ST?? 

|I|<4 

$15 4 < |I| < 8 

f :i 8 < I 1 ! 

;j It is noted that the references above to MSB, LSB and MID. SIG. B are to the most, least 

{ IJ and middle significant bits, which are each pairs of bits that make up the constellation. 
T%20 Further, the Q values are the same as the I values provided above. 

; 1 These metrics have a number of advantages. First, they are able to provide soft decision 

values to data that may not be able to be directly decoded via soft decision decoding using 
Euclidean distances, providing an advantage over the alternative of hard-decision decoding. 
Secondly, these metrics provide a more extreme weight when the position is further from the 

25 transition region. For example, in the absence of noise, a 64-QAM I-channel input of 1, 3, 5, 
or 7 corresponds to a 1 for the msb. However, an input of 1 is on the transition edge, and 
could be noisy, so it is given a lower weight (0.75) compared to anything greater than a 2, 
which is far enough away from the transition edge to be more certain of the value of the msb. 
This allows the algorithm to provide more weight in decoding where there is more certainty, 

30 providing a greater chance of correctly decoding all of the bits. 

These metrics can be weighted further based on knowledge of the probability that the 
symbol is in the presence or absence of noise or interference of any kind. This degradation 
may be due to any effect in time, space, or frequency, such as an interferer or a poor channel 
response in a given frequency bin for an OFDM signal, or an intermittent interferer arriving 



LSB 
|I| * 0.25 
(8 - |I|) * .25 
0 
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in given symbols in time, or a given antenna receiving a weak signal in space, or any 
combination of any effect in these three independent domains. For example, in the 802.1 la 
instantiation of this discussion, we may use our channel estimate to detect frequency bins 
which may possess low SNR due to interference or channel response. All soft decisions from 
5 these frequency bins can then be weighted according to a predetermined estimate of their 
reliability, giving more weight to the reliable bits and further increasing the ability to decode 
all bits correctly. 

The bit values calculated above are then used in the Viterbi decoding algorithm similarly 
to the above — the bit values are used as received binary subsequences, and metrics of trellis 
10 branches and states are calculated therefrom. Consider the example of FIGs. 2 and 3, where 
in place of C(X) = (1 11010 00 010010101 1) we have a 16QAM system in which the 
received (I, Q) pairs are [(-0.5, +1.6) (-3.4, -1.2) (-2.7, +3.2) (-0.7, -1.4) (-2.5, +0.9)], which 
;'P generally corresponds to the same sequence with some noise added. Applying the above 
i transfer functions, the following "bit" sequences are obtained: 
ll5 (0.7, 1,1, 0.25) corresponding to 0 1 1 1; 

f (0.9, 0, 0, 0) corresponding to 0 0 0 1 ; 

r (0, 1 , 0. 1 5, 0) corresponding to 0 0 1 0; 

*J (0.8, 0, 1 , 0. 1 5) corresponding to 0 1 0 1 ; and 

;ij (1, 0.95, 0.25, 0) corresponding to 0 0 1 1. 

; 520 Applying these figures to the trellis diagram of FIG. 2C obtains the weighted trellis 

* shown in FIG. 9. By working backward from the least-cost path in the last stage as above, 

the best guess state transition path of the coder 10 can be obtained, and from that the original 

input sequence as before. 

It should be noted that although the above computations have been illustrated using the 
25 decimal system for ease of illustration, in practice it will likely be more convenient to 

calculate the partial bit values as n-bit binary values, with each binary step corresponding to 

an increase of l/2 n from the next lower value. In the preferred embodiment n=4~(5), which 

was chosen as a tradeoff between hardware complexity and decoder performance; other 

values will also prove workable. 
30 So, a bit value of 0.5 would be represented by a value of 10000, a value of 0.8 might be 

rounded up to 1 1001 (0.8006), etc. Slightly larger accumulators can be used to store the state 

metrics; a six bit accumulator could be used to store any of the state metrics shown in FIG. 9. 

This will simplify circuitry implementing the embodiment immensely and speed it up a great 

deal. 
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Generally speaking, more than six bits might be required in the accumulator, especially 
when more than nine bits are decoded. Also, those of ordinary skill in the art will recognize 
that techniques are known for limiting the required number of state metric bits, even for 
frames of arbitrary length. 

It will also be appreciated by those skilled in the art that a number of different 
fabrication technologies may be used to implement the preferred embodiment — lower 
throughput systems may be implemented in software, while higher throughput systems may 
use a hardware-oriented approach such as a custom chip, an FPGA, a gate array or the like. 

A second preferred embodiment of the present invention is similar to the first in its use 
of many- valued bit values; however, in the second embodiment the transitional portion of the 
curves are not linear but curved. FIG. 10A shows the transfer curve for bo and b2, and for bj 
and b3, corresponding to FIGs. 7A/8A and 7B/8B, respectively. The transitional portions of 
these curves generally can be expressed either in formulae such as those below (using the Q 
parameter for exemplary purposes) or, in practice, in a lookup table. 

For FIG. 10A, 
Q < -1 -» bit = 0 

1 O 2 

2 2 

1 o 2 

2 2 

Q> 1 -+bit = l 
For FIG. 10B, 

Q<-3 -> bit = 0 

-3< G ^-2-,^=I-fe±2)l 
* 2 2 

-2^<-l-^to=I + M±2)l 
2 2 

-1 < Q < 1 -» bit = 1 

2 2 
* 2 2 
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Q > 3 bit = 0 

Of course, these formulae are only exemplary, and a wide variety of non-linear transi- 
tional curves may be used in their place. For example, having a steeper transition near the 
5 flat parts of the functions can produce slightly better results than a linear transition region. 

Generally, an embodiment according the invention described above allows a receiver to 
independently differentiate predictably reliable data from predictably unreliable data. Both 
single carrier and multicarrier bit transmissions often provide sporadic, unreliable data. For 
example, in a single carrier system in which the channel bandwidth is narrow within the 
1 0 context of a selected bit transmission rate, it is often difficult to accurately receive sequences 
containing long strings of uninterrupted 1 ? s or 0 ? s. Similarly, bits transmitted through fading 
channels in a multi-carrier system can be problematic, as can bits contaminated with frequen- 
ce cy dependent noise (i.e., bits transmitted at lower frequencies are less reliable than those 
f i transmitted at higher frequencies) and subject to other non-ideal transmission conditions. 
;U 5 FIG. 1 1 illustrates an example of the distinction between good and bad data in a multicarrier 
i iA environment. 

s The following is a discussion of the preferred puncturing method and system within the 

context of a multi-carrier environment that allows for the differentiation noted above. 
t= Assume a Coded Orthogonal Frequency Division Multiplexed ("C-OFDM") transmission 

."^20 within a wireless LAN environment, such as the environment provided by the IEEE 802. 11a 
protocol. In such a C-OFDM environment, such as illustrated in the transmitter 100 shown in 
Figure 12, blocks of input bits are first coded using a coder 120 (redundancy is added to the 
bits to correct future errors) and then interleaved with interleaver 130 (the bit order is 
changed in a known, pseudo-random manner). The coded and interleaved block is then 
25 passed to a modulator 140, which maps the block to a constellation point using, for example, 
gray codes. Once modulated, the modulated signal is passed to an Inverted Fast Fourier 
Transform ("IFFT") module 150 as illustrated in FIG. 12. Every nth bits from the block 
(such as 1, 2, 4, and 6 bits for BPSK, QPSK, 16QAM, and 64QAM modulation techniques, 
respectively) modulates one of forty-eight available carriers. The multicarrier symbol is then 
30 modulated to the RF frequency using a conventional D/A and RF amplification circuit 1 60 
and transmitted over an antenna and received at the receiver 200, which is illustrated in more 
detail in FIG. 13, described below. 
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For the sake of the provided example, a X A rate coded stream will be assumed, wherein 
the coder adds a redundant bit for every input bit, thereby doubling the number of bits 
transmitted for a given amount of data. In order to reduce part of this overhead, some of the 
coded bits may not be transmitted - that is, they may be hard punctured at the transmitter - 
5 providing, for example, a resulting % rate code, as is conventional. However, the Viterbi 
forward error correction ("FEC") decoder, unaware of any hard puncturing that may have 
occurred on the transmitter side, decodes the received bits as if they were transmitted as Vz 
rate code. The hard punctured bits are therefore added to the received bits at the receiver 
before being passed into the Viterbi decoder and decoded as normal Vz rate code. However, 
10 the Viterbi decoder is also given the location of the hard punctured bits so that those bits are 
not used in calculating the norm of the difference between the received bit and the expected 
bit on the Viterbi trace (generally known as the distance metrics). In this case, all distance 
v!) metrics are zeroed out (hard punctured). The hard punctured bits will therefore have no 
1 J actual effect on the decoding process, and the performance of the error correction capability 
;; 15 of the Viterbi depends only on the non-punctured bits. This hard puncturing operation is 

performed at the transmitter to allow a higher data rate of transmission at the expense of less 
redundancy. 

! * The present invention can implement at the receiver the soft or hard puncturing of bits 

that were transmitted but not deemed as worthwhile as other bits using the same mechanism 
::20 for the bit metric puncturing or weighting that is used when hard puncturing from the 
•i i * transmitter, as described herein. Thus, in a preferred implementation of the present 

invention, any bit can be flagged to be soft or hard punctured. The Viterbi decoder can 
discard flagged bits from its metric calculations (hard puncturing) or de-emphasize certain 
bits by weighting the associated distance metric (soft puncturing). Through these methods, 
25 better bits get higher weighting in the Viterbi decoding process (distance metric calculations), 
and worse bits either get lower weights or get punctured (if below a threshold, which could 
be dynamic). 

A further description of the receiver 200 that can provide for both soft and hard 
puncturing according to the present invention will now be described with reference to Figure 
30 13, with particular attention provided to those modifications necessary to a conventional 
receiver to implement the present invention being described. Hard puncturing will not, 
however, be described differently than soft puncturing, since hard puncturing as used herein 
is a degenerate case of soft puncturing, in which a weight is set to zero and hard puncturing 
of the bit metric results. 
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As shown, the receiver 200 includes a front end/analog to digital converter 212 and 
timing recovery circuit 214 and fast Fourier transform (FFT) circuit 220 that operate in a 
conventional manner. Further, as described herein, the channel correction circuit 222, 
deinterleaver circuit 224 and Viterbi decoder 226 also operate in a conventional manner, 
5 except that the inputs received by these circuits are modified due to the establishment of the 
weighting values described herein that are used to implement the soft puncturing that takes 
place. 

An overall description will initially be provided. As shown, the output from FFT 220 is 
also input into a dynamic channel weighting circuit 230. The dynamic channel weighting 
10 circuit 230 provides a dynamic weighting, as described herein, that is associated with 

dynamic characteristics based upon dynamic changes in channel performance, for each of the 
bins in the channel. This dynamic weighting is combined with a static weighting, obtained 
v [!: from a static channel weighting circuit 250, for each of the bins in the channel, in a combiner 

260. The resulting channel weight is then used to alter the value of the bit metrics used by 
:^,15 the Viterbi decoder 226 when making branch metric calculations, 
i** Further details will now be provided. The dynamic channel weighting circuit 230 

fg includes a channel estimator 232, which, as known, can provide an estimate of the 

;! performance of each of the bins in the channel. 

j The channel estimate for each of the bins is output to channel inverter 234, which 

^20 obtains an inverse of the channel estimate, which inverse can then be applied to the channel 
* correction circuit 222 in order to assist in improving the signal. 

The output from the channel estimator 232 is also applied to a normalization circuit 236. 
The normalization circuit normalizes the signal for each of the channel bins in some manner, 
to determine whether the signal has a sufficient strength that suggests it should be fully 
25 weighted by the Viterbi decoder, or has some lesser strength, which thus warrants a dynamic 
weighting value that attributes less significance to the particular signal. The channel estimate 
is used to determine the strength of the signal, which is then compared to the normalized 
signal. The signal strength is typically obtained in terms of the power or magnitude of the 
signal corresponding to each bin in the channel, as output from the FFT 226. The usage of 
30 power, calculated as I 2 + Q 2 (which can be derived from the I and Q obtained as a result of 
the channel estimate) for the complex channel estimate for each bin, is preferred, since to 
calculate a magnitude, a square root function is required to calculate an absolute value, 
thereby simplifying the calculation, since there is a 1 to 1 mapping of power to magnitude. 

! 
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A specific example of dynamic weighting values for various thresholds is provided in the 
table below: 



Weight 


Magnitude Threshold 


Power Threshold 


1 


V2/2 


Vi 


V2 


V2/4 


1/8 


% 


V2/8 


1/32 


1/8 


V2/16 


1/128 


1/16 


V2/32 


1/512 


0 


N/A 


N/A 



For proper interpretation, any bin with a power estimate between 1/4 and Vi of the max 
5 magnitude gets weighted by V2 in the Viterbi decoder 226, assuming that there is no static 
weighting component, as described hereinafter. 

Thus, in a preferred embodiment of the invention, different dynamic weighting values, 
'tt which can be represented by three bits, are used. The dynamic weighting values that are less 
: A than 1 are derived based upon an amount that the signal deviates from some normalized 
r ilO signal value. This normalized signal value could be, for example, the maximum power or 
j f magnitude or the mean power or magnitude of the signals in the channel. In any event, for 

each signal value that is further from the normalized signal value, that signal value is 
jrl associated with a progressively lesser weight, as illustrated above. Accordingly, based upon 
uj the strength of the signal, a corresponding dynamic weighting value is derived, which is then 
; 5 1 5 stored in a channel memory 238. 

For each of the bins in the channel, the corresponding dynamic weighting value is then 
output from the channel memory 238 and input into the combiner 260. 

The combiner 260 will combine, in a manner more fully described hereinafter, each 
dynamic weighting value with a corresponding static weighting value obtained from static 
20 channel weighting circuit 250. 

Static channel weighting circuit 250 will now be described in further detail. Before 
providing such further description, however, it is noted that the term static is used with 
respect to this weighting circuit not to suggest that these weights do not change, but to 
distinguish these weights from those that are obtained through usage of the channel estimates, 
25 as described above. Further, typically, the nature of the interference which results in the 
static channel weighting values will occur for a longer period of time, though that does not 
necessarily need to be the case. An example of an interferer that may cause a static weight 
that will lower the significance of the signal in the Viterbi decoder 226 is circuit noise at a 
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predetermined frequency corresponding to a particular bin within the channel currently being 
used. 

As illustrated, static channel weighting circuit 250 includes a static weighting memory 
252 that stores a channel performance table. In the preferred embodiment, the channel 
5 performance table includes initial static weight values for each bin of each channel of 
interest. Each initial static weight value for each bin can be influenced by various factors, 
some of which may be entirely static, but many of which will change over time. Examples of 
factors that can change these initial static weight values are discussed further hereinafter. 
For a given bin in a given channel that has its dynamic weight value read out of channel 
10 memory 238, a corresponding initial static weight value will be output from the memory 252 
to interpolator 254. Interpolator 254 will adjust the initial static weight value for specific 
^ predefined operational characteristics. For example, the interpolator 254 may be provided 
P with the rate at which data is being transmitted, the determined frequency offset, as well as 

the analog front end gain currently being used. Interpolator 254 will then adjust the 
rLl5 corresponding initial static weight value to determine an actual static weight value. The rate 
! ^ and analog gain settings used will either turn on or off the soft puncturing of certain bins, or 
will scale the amount of puncturing. The estimated frequency offset will be used to shift 
which bins are punctured. Interpolator 254 will also preferably scale by a power of 2 when 
: y obtaining the actual static weight value. Thus, while the interpolation is somewhat coarse, 
T20 with a thresholded frequency offset, and a predefined, power of 2-based interpolation 
= r « function, a more complex scheme is not necessary. 

As noted above, the actual static weight value is then combined with the corresponding 
dynamic weight value, for each bin of the channel. In the preferred embodiment, with the 
specific weights discussed above, the combiner 260 is implemented as a look-up table 
25 emulating a simple multiplier. The resulting bin weight value is then output to the 

deinterleaver, which will store the resulting bin weight value along with the corresponding 
deinterleaved data signal in a buffer contained therein, so that both can be output at the same 
time for usage by the Viterbi decoder 226, as discussed generally above and which will be 
described more fully below. The deinterleaved data signal contains the I or Q value, as 
30 appropriate, for a given coded bit, as well as the bit position (msb, lsb) of that bit in the 
constellation to determine which bit metric table to use. 

Within the Viterbi decoder 226, a portion of which is illustrated in Figure 14, the 
deinterleaved data signal is input into a bit metric table 300, which, for the signal that has 
been input, determines a bit metric estimate value, which is an estimated probability that a 
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particular bit is a "1 A complement bit metric estimate value, which is an estimated 
probability that a particular bit is a "0," is also obtained using a l's complement circuit 302. 

In a preferred embodiment of the invention, each bit metric estimate value output of the 
bit metric table will be a five-bit number, such that the probability can be characterized as 
5 one of 32 different values. Each bit metric value and corresponding complement bit metric 
value is then operated upon using the corresponding resulting bin weight value that has been 
established. Accordingly, the bit metric values are modified by combiners 304 and 306 using 
the weight value, also referred to as being soft-punctured, preferably by shifting the bit metric 
values by an amount determined by the weight value. Thus, in subsequent branch metric 
10 calculations, the weights used by the Viterbi decoder 226 take into account both the channel 
and other a-priori characteristics that allow differentiation between good and bad bits. 
Within the Viterbi decoder 226, the bit metric values, or decision metrics, are used in 
w determining how to proceed to the next state of the trellis, with each of the decision metrics 
^ being modified by the puncturing that has occurred. Thus, if hard puncturing has occurred, 
j 5 thereby giving no weight to the punctured bit, each state in the trellis will simply select the 
^ previous possible state with the lowest state metric and maintain that metric. If soft 
] * puncturing has occurred, thereby giving some weight to the punctured bit, each state in the 
* trellis will simply select the previous possible state taking into account the weighted decision 
: ;j metric. 

1;20 The above description thus provides a description of how to implement soft and hard 

puncturing in the receiver according to the present invention. The following description will 
now provide a further description of the manner in which the soft and hard puncturing in the 
receiver, according to the present invention, are useful. 

Thus, as one example, puncturing can be used to avoid known jamming signals. Thus, 

25 as described above, if particular frequency bins in an OFDM based system (or other samples 
in a single carrier system) are thought to be erroneous, they can be punctured to be ignored 
completely, or weighted less to reduce the negative effects of errors while allowing some 
useful information to be extracted. Particular frequency bins could be erroneous due to 
transmitter or receiver circuit noise, jammers occupying the same frequency space from 

30 incompatible communication protocols, or devices emitting RF interference such as 
consumer microwave ovens. 

To address frequency-hopping interference (for example, Bluetooth signals, wherein the 
interference appears pseudo-randomly in multiple frequency bins), one could either note 
which frequencies are being trampled by detecting low-SNR, or else provide just enough of 



i 
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the necessary jammer radio circuitry to determine the hopping sequence. For example, to 
mitigate 2 A GHz Bluetooth interference in a high-rate 2.4GHz IEEE 802.1 lb system, one 
could build enough of the Bluetooth radio and protocol stack to lock on to hopping 
sequences. This hopping-sequence lock would then allow a prediction of which bins to 
puncture (completely ignore) or weight less (partially ignore). 

Another method to detect jammers or interferers in an OFDM system with channel 
estimates (such as required in IEEE 802.1 la) is to examine the channel estimates and detect 
abnormally high or low carrier bins with respect to their neighboring carrier bins. For 
example, if three bins show a much greater channel response, particularly with a greatly 
different phase, it is likely that they are from a jammer. This is because there is a certain 
degree of smoothness anticipated in the channel response. If two carriers reveal phase 
differences in the two long training symbols, this too indicates that they may not be from an 
expected protocol-compliant signal, in which case such signals can be punctured. 

Also, maximum likelihood (ML) criteria for soft puncturing may be used. The optimum 
weighting based on ML criteria is proportional to the inverse of the channel estimate. 
Therefore, for soft puncturing, we use the inverse of the channel estimate to weight the 
Viterbi metrics, providing a metric to quantify the reliability of received bits, which, in turn, 
correspond to frequency bins. In a preferred embodiment, the invention includes the ML 
criteria to quantify frequency bin reliability. The first step in the process is to define the 
likelihood function L(X | 0). For example, for a noise process that is cyclo-stationary and 
Gaussian, the likelihood function at the receiver can be written as; 

L{X\&) = Ke a ' 2 (1) 
where K is a constant and i represents the frequency bin. The likelihood function shows how 
likely it is to receive a symbol X given that symbol 0 was transmitted, where C\ is 
proportional to the inverse of the magnitude of the channel estimate. 

The log of the likelihood function, which provides an estimate of the likelihood of 
receiving a particular sequence where Gaussian-distributed noise is found, would be: 

L(X 1 0) = -log(jQ£ (X ' ~P' )2 (2) 

The Viterbi distance metric corresponds to the log of the likelihood function. From the 
above it is apparent that the distance metric should be proportional to <7i, where <J\ is a 
function of the channel estimate. It should be noted that the argument above was based on an 
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"equivalent noise" analysis, where it is assumed the received frequency domain data is scaled 
by the inverse of the channel. 

The frequency domain system can therefore be written as 

Yi^HiXi + Ni (3) 
5 where Hi is the channel at bin i. Now consider Hi as the channel estimate. Scaling by the 
inverse of the channel estimate leads to 

Y i = (Hr 1 )Y i ^X i + (H i - 1 )N i (4) 
where Oi 2 = |Hi _1 |a 2 . This shows the amount of noise enhancement in each frequency bin i , 
an amount proportional to the channel inverse. 
10 Thus, according to the formulated likelihood criteria, the Viterbi distance metrics are 

preferably scaled by the above amount for a given SNR estimate. So for a channel H(k) 
where k represents the frequency bin number, 
5 1 

f= new branch metric (k) = branch metric (k) • — (5) 

£ K 

r ; t max I H k I 

: 3 where a t = - 



1 5 For comparison to the background art, note that the known log likelihood function for a 

coded OFDM system with channel Hjk (j=time, k=bin) as described in Alard et ah, "A New 
System of Sound Broadcasting to Mobile Receivers" (1988) is 

Y\ Y jk-~ H jk* C jk 
2cT 



MY | C) = *" ^ 2 (6) 



For BPSK and QPSK, the available art describes a decoding method. For higher order 
20 modulations, such as 16-QAM, 64-QAM and others, the "metric" term 

M jk =\Y jk -H Jk *C Jk ( (7) 

is typically separated into per-bit metrics that approximate the L2 Euclidean metric but is 
easier to calculate. This means that the "sample" or "symbol" metric is composed of a sum of 
bit metrics, e.g. 

25 M Jk =^m[lW jk M Jk ,C Jk ) (8) 

where m is some specific per-bit metric function. In a Gray-coded 16-QAM system a 
particular choice of per-bit metric functions is described in "OFDM for Wireless Multimedia 
Communications", van Nee and Prasad (2000). 
Equation 6 can also be written in the form: 
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I 



2G[j,k? 



(9) 



l Jk " J 

where Hjk can be replaced with an estimate H jk . For implementation purposes, this is a more 
amenable form that, for example, permits the use a scale-independent lookup-table for the 
pre-bit metrics. 

5 The inventive method comprises of extending the available method to: 

1. Use the formulation of Equation 9. 

2. Compute the numerators of Equation 8 using any type of separable, per-bit metric. 

3. Apply an estimate H jk of the channel H jk to "weigh" the denominators as indicated in 
IT eqn4. 

( : ?10 This is an extension of the BPSK and QPSK scheme described in the available art to 

!7 higher order modulation configurations that have multiple bits per I and Q channel, 
w In an OFDM signal, stationary background noise and a stationary or slowly varying 

[ channel gives rise to a situation where Hjk is periodic in the bin number k (same value in 

same bin for consecutive symbols), so Hjk=Hk. This has the same effect as n-equivalent 
. J 5 cyclo-stationary noise 

<T'(k) 2 = ° 2 , (10) 

- \\H k f 

along with a flat equivalent channel model. 
An additional optional feature of the inventive method and system entails using 
decoded/error-corrected parts of a data packet to train or improve decoder performance for 
20 the rest of the packet. Thus "decisions" about earlier parts of a packet "direct" the decoding 
of later parts of the packet. Used iteratively, decisions made on an initial pass direct the 
decoding of later passes. 

Of course, an iterative implementation can lead to pitfalls. For example, an error 
occurring early in the packet or sequence causes false training later in the packet or sequence, 
25 thus worsening performance. This is known generally as error propagation. However, in the 
case of protocols such as 802.1 la, if there are any errors in the packet, the entire packet is 
typically discarded. Thus, all post-Viterbi decisions are essentially "correct" with respect to 
training for later decisions. While the iterative approach can also be applied to pre-Viterbi 
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decoder data (to reduce latency), care must be taken to avoid excessive error-propagation, 
since pre-Viterbi data can be incorrect yet still result in a correct packet. 

Decision directed feedback in an OFDM-based receiver, such as used in 802.1 la- 
compliant systems, can operate by passing the output of the Viterbi decoder through the 
5 transmit chain up to but not including the IFFT. Thus, the decoded data is fed through the 
FEC coder, interleaver and mapper and represents the data that is expected in the receive 
chain after the FFT and channel correction. There will be a latency between when the data 
enters the receiver and when it exits the Viterbi decoder due to the processing delay of the 
demapping, FEC decoder, FEC encoder, and interleaver. The post-FFT receive data must be 
10 delayed similarly so that it will be available for comparison at the same time that the 
reconstructed data is ready. FIG. 15A illustrates the described apparatus. 

There are cases where decision feedback prior to decoding may be preferable. This may 
; : : be the case if the latency of waiting for the decoding and re-encoding is harmful, such as in 
"Q the case of short packets, or if the added circuitry to perform re-encoding is overly costly. In 
;[15 this case, if the occasional error propagation described above has a tolerable effect, decision 
r * feedback prior to decoding may provide a better solution to estimate the per-bin SNR to the 
™ Viterbi decoder weighting circuitry. FIG. 15B illustrates the described apparatus. 
O Once the reconstruction has been performed, a stream of received and desired frequency 

; y domain data is made available that can be used to adapt the receiver. Essentially, all of the 
;20 subcarriers are made to act like pilots, as their desired and actual data is known. This is useful 
in determining per-bin SNR for use in selective puncturing/weighting, in which the signal-to- 
noise ratio per sub-carrier (or FFT bin) can be determined by computing the standard 
deviation of the desired frequency domain data and dividing it by the standard deviation of 
the difference between the desired data and the received data. The initial SNR estimate based 
25 on the channel estimate magnitude can then be adapted or replaced with the new normalized 
SNR estimates from the decision feedback circuitry. Selective puncturing or weighting of 
bins can then be executed in a manner based on the determined SNR, as described earlier. 
Bins with smaller SNR would be punctured or weighted less than bins with greater SNR. 
The SNR estimates can be used on a per-packet basis and/or carried over between 
30 packets. This carrying over may be gated by the determined success or failure of the packet 
based on a checksum. Carrying over SNR estimates allows for greater averaging and is more 
effective acting for earlier parts of the packet, but care must be taken since two packets can 
come from different sources, under different channels, and at different times. Thus, some 
estimates from earlier packets will not be valid for later packets. Thus reuse of data between 
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packets is most useful for receiver-specific processing, such as spurious tone estimation and 
long-term jamming detection and estimation. Jammer estimation may or may not vary 
between packets. 

The present invention has been described above in connection with a preferred 
embodiment thereof; however, this has been done for purposes of illustration only, and the 
invention is not so limited. Indeed, variations of the invention will be readily apparent to 
those skilled in the art and also fall within the scope of the invention. 
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